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Control over the quantum dynamics of chaotic kicked rotor systems is demonstrated. Specifically, 
control over a number of quantum coherent phenomena is achieved by a simple modification of the 
kicking field. These include the enhancement of the dynamical localization length, the introduction 
of classical anomalous diffusion assisted control for systems far from the semiclassical regime, and 
the observation of a variety of strongly nonexponential lineshapes for dynamical localization. The 
results provide excellent examples of controlled quantum dynamics in a system that is classically 
chaotic and offer new opportunities to explore quantum fluctuations and correlations in quantum 
chaos. 
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I. INTRODUCTION 

The quantum kicked rotor (KR) and its classical ana- 
log, the standard map, have long served as a paradigm 
for quantum and classical chaos pj. Due to its atom 
optics realization 0, the KR has recently attracted re- 
newed interest. The KR is also of considerable interest in 
a variety of other fields such as condensed matter physics 
0, , molecular physics 0, |(J , and quantum information 
science 

One well-known quantum effect in KR is "dynami- 
cal localization" (DL) 1]. That is, although a clas- 
sical kicked rotor displays unrestricted diffusive energy 
increase due to classical chaos, only a finite number of 
rotational states are excited in the quantum dynamics, 
with the quantum excitation probability versus the rota- 
tional quantum number typically assuming a character- 
istic exponential lineshape. DL is a pure quantum coher- 
ence effect and is therefore very sensitive to decoherence. 
For example, previous studies have demonstrated that 
noise |0j, lift llf| , nonperiodicity in the kicking sequences 
[l2Lll3| |. and quantum measurements Q can destroy DL. 

As a coherence effect, DL is also indicative of the possi- 
bility of controlling the KR dynamics via quantum effects 
0,0)0. Indeed, we recently showed that the quan- 
tum phases describing the initial rotor quantum state 
can be manipulated to effectively control quantum fluc- 
tuations in quantum chaos and thus enhance or suppress 
quantum chaotic diffusion [l7L llgj in KR. However, ma- 
nipulating quantum phases in initial states cannot change 
the unitary evolution operator of the system. Hence, nei- 
ther the average dynamical localization length nor the 
characteristic lineshape of dynamical localization can be 
altered in this way. 

Motivated by interest in controlled classically chaotic 
quantum dynamics |17lll8lll9| . and to gain more insights 
into the nature of DL, we consider controlling DL and the 
associated energy absorption via a modified kicked rotor 
(MKR) system, in which the phase of the kicking field, 
or the timing of the kicking sequences (hence the evolu- 
tion operator), is actively manipulated. In particular, we 



consider an MKR system in which the sign of the kicking 
potential is periodically changed, or alternatively time 
delayed, after a certain number of kicks. As shown below, 
such a slight modification of KR has profound effects on 
the dynamics: whereas periodically changing the sign of 
the kicking potential does not destroy DL, it dramatically 
changes the quantum diffusion dynamics of KR as well as 
the nature of DL, offering entirely new opportunities for 
controlling dynamics, as well as understanding quantum 
fluctuations and correlations in quantum chaos in peri- 
odically kicked systems. For example, we demonstrate 
that, with the sign of the kicking potential of KR peri- 
odically changed: (i) the dynamical localization length 
is significantly increased so that the energy absorption 
is strongly enhanced; (ii) classical anomalous diffusion 
(which can be slower than quantum anomalous diffusion 
under certain conditions pd|) can enhance control even 
when the effective Planck constant is about an order of 
magnitude larger than the relevant classical phase space 
structures, and (iii) DL may display strong deviations 
from purely exponential lineshapes. 

This paper is organized as follows: In Sec. [H]we intro- 
duce the modified kicked-rotor model. We then present 
results in Sec. HHI on the enhancement of dynamical local- 
ization length, with qualitative explanations based upon 
a known result from the Band Random Matrix theory 
[Ulilli^. In Sec. |TV|we show control of DL in a dif- 
ferent regime, where the dynamics can be tied to a new 
mechanism, i.e., the creation of new structures in classical 
phase space. Section |3 contains the results on coherent 
manipulations of the lineshapes for DL. We conclude the 
paper with a brief discussion in Sec. IVII 



II. A MODIFIED KICKED-ROTOR MODEL 

The KR Hamiltonian is given by 
H KR (L, 6, t) = L 2 /2I + A cos(0) ^ 5{t/T - n), (I) 
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where L is the angular momentum operator, 9 is the con- 
jugate angle, I is the moment of inertia, A is the strength 
of the kicking field, and T is the time interval between 
kicks. The basis states of their Hilbert spaces are given 
by | to), with L\m) — mh\m). The quantum KR map 
operator for propagating from time (N — + )T to time 
{N + 1 - 0+)T is given by 



(MKR) system whose Hamiltonian is given by: 



F K r = exp 
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exp[— ifccos(0)], 



(2) 



with dimensionless parameters k — XT/% and the ef- 
fective Planck constant r = KT/I. For later use we 
also define the dimensionless scaled rotational energy as 
E = (L 2 )t 2 /2Ti 2 , where (■) represents the average over 
the quantum ensemble. In the \m) representation, Fkr 
takes the following form [23| : 

( mi \F KR \m 2 ) = exp (ilm?) t mi ~" 12 J mi _ ma (fc), (3) 

where J mi - m2 (k) is the Bessel function of the first kind 
of order (mi — 772,2). 

The classical limit of the KR quantum map, i.e., the 
standard map, depends only on one parameter k = kr 
and is given by: 



L N = L N -i + Ksin(0jv-i); 
On = ®n-i + Ln, 



(4) 



where L = Lt/K is the scaled c- number angular momen- 
tum and (£jv, &n) represents the phase space location of 
a classical trajectory at (N + 1 — + )T. For later discus- 
sion we note that for particular values of k, the classical 
map Eq. Q can generate accelerating trajectories whose 
momentum increases (or decreases) linearly with time (at 
least on the average). These trajectories are called trans- 
porting trajectories [24|. To see this consider the initial 
conditions: (L = 2ttIx,6 — ±vr/2) for k = 2nl2, where 
li and I2 are integers. Clearly, these phase space points 
are shifted by a constant value (±2^/2) in L after each 
iteration, resulting in a quadratic increase of rotational 
energy. These transporting trajectories are rather sta- 
ble insofar as they may persist for values of k close to 
27r?2 (with their average momentum shift after each iter- 
ation oscillating around the constant value ±271^2), thus 
giving rise to transporting regular islands |24| . i.e., the 
accelerator modes in the KR case. If classical trajecto- 
ries are launched from the accelerator modes, they simply 
jump to other similar islands located in adjacent phase 
space cells. For trajectories initially outside the accelera- 
tor modes, the "stickiness" of the boundary between the 
accelerator modes and the chaotic sea induces anomalous 
diffusion over the energy space, i.e., energy increases in a 
nonlinear fashion, but not quadratically. This is intrinsi- 
cally different from the case of normal chaotic diffusion in 
which energy increases linearly with the number of kicks. 
We introduce here a slightly modified kicked rotor 
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(L, 0, t) = L 2 /2I + A cos(0) fM(n)6(t/T - n),(5) 



where /jf(f») is real, |/m(?i)| = 1, and Jm(ji) changes 
sign after every M kicks. That is, the only difference 
between KR and MKR is that in the MKR the sign of 
the kicking potential is changed after every M kicks. 
The effect of changing the sign of the kicking poten- 
tial can be further understood in terms of the time- 
evolving wave function, which can be expanded as a su- 
perposition of |m) states: ^2 m C m {0\m) , with the ex- 
pansion coefficients C m . Since cos((9 + w) = — cos(0), 
and Y.C m {0 + K\m)=J2 m (~l) m C m {0\m), changing the 
sign of the kicking potential is seen to be equivalent to 
adding a n phase difference between all neighboring ba- 
sis sates. Compare this now to the effect of a time de- 
lay td — 2-kT/t between two neighboring kicks. Due 
to the free evolution of the rotor any two angular mo- 
mentum eigenstates |m) and \m + 1) will acquire in 
time td an additional relative quantum phase given by 
exp{iTt d [(m + I) 2 - m 2 ]/(2T)} = exp[i(2m + 1)tt] = 
exp(i7r). As such, the MKR can be also realized by in- 
troducing the time delay td after every M kicks. 

For times (N - 0+)T to (N + M - 0+)T, the MKR 
quantum propagator can be written as 



F, 



MKR 



d 2 

° X 1' \ l7r Q02 ) F KR 



dfMr- 



(6) 



where this equation defines D as the free evolution op- 
erator over time td- Here F^ R denotes M applications 
of Fkr- From Eq. ©, one sees that the only differ- 
ence in time propagation between KR and MKR for ev- 
ery M kicks is the D operator, whose matrix elements 
(mi I D I m.2) are given by 



(mi|D|m 2 ) - (~l) mi S mim2 . 
Thus we have 

(miliWilK) = (-l) mi (m 1 |F^|m 2 ) 



(7) 



(8) 



Note also that the classical limit of the MKR quantum 
map [Eq. ©] is given by 



L N+ i = L N + Kf M (N)sm(6 N ), 
On = On + Ljf+i. 



(9) 



With respect to model systems in the literature, the 
MKR here can be regarded as a specific realization of 
the so-called generalized kicked-rotor model, which was 
first introduced in Ref. 25], in the context of quantum 
anti-resonance. However, it is very different from the 
amplitude-modulated kicked rotor systems [12l |26( pre- 
viously studied because (i) the kicking field strength here 
remains constant, and therefore any interesting results 
arise from pure phase modulations; and (ii) as shown be- 
low, the dynamical localization is significantly altered, 
but not destroyed. 



3 



III. ENHANCED DYNAMICAL 
LOCALIZATION LENGTH 

As an example, consider the case of k — 4.0, r = 2.0, 
and M = 50, whose classical limit for both the KR and 
MKR is fully chaotic and displays normal chaotic diffu- 
sion. We demonstrate below that the dynamical localiza- 
tion (DL) and the related energy absorption associated 
with Fmkr are dramatically enhanced over that associ- 
ated with Fjf R . 
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FIG. 1: Phase control of dynamical localization achieved by 
changing the sign of the kicking potential after every 50 kicks. 
k = 4.0, r = 2.0, and the initial state is |0). (a) A comparison 
between KR (the narrow lineshape, solid line) and MKR (the 
broad lineshape, dashed line) in terms of the probability P(m) 
of finding the system in the state \m) after 4 x 10 kicks, (b) 
The time dependence of the dimensionless scaled rotational 
energy E in each case of KR (solid line) and MKR (dashed 
line). Note that the solid line lies very closely to the E — 
axis. 

Figure displays the angular momentum distribu- 
tion P(m) after 4 x 10 5 kicks, starting with the initial 
state |0), for both KR and the MKR with M = 50. 
The exponential line shape of P(m) shown in Fig. {T^l 
indicates that DL occurs in both cases. By fitting 
P(m) with exponentials P(m) ~ exp[— \m\/lj(ii] ancl 
P(m) ~ exp[— \m\/lMKR\ for KR and MKR, respec- 
tively, one obtains that the dynamical localization length 



Imkr ~ 140.0 is significantly larger than Ikr ~ 7.0. This 
clear difference in dynamical localization length is also 
reflected in the energy absorption shown in Fig. hi 
particular, while the energy absorption of KR saturates 
after a few kicks, the MKR system continues to absorb 
energy in a more or less linear manner for as long as 10 4 
kicks. Enhancement of dynamical localization length and 
energy absorption is also observed for other values of M, 
and for a wide range of parameters k and r. 

As is well known, the DL of KR can be traced back to 
the localization properties of the eigenstates of the quan- 
tum map operator [Eq. J2J]. Note first that, due to the 
rapid decay of J mi - m2 (k) with increasing \mx — and 
the pseudo-random nature of the function exp(irr7i 2 /2) in 
mi, the quantum map operator Fkr in general assumes 
a band structure and behaves in a pseudo-random man- 
ner in the |m) representation. Hence, below we qualita- 
tively consider the MKR results in terms of a well-known 
feature from Band- Random-Matrix theory 0, 1221 123| , 
namely, that the larger the band width of the quan- 
tum map operator the larger the dynamical localization 
length. 

Note first that the matrix (m\FKR\m') is pseudo- 
random. However, we do not expect {m\F^ R \m') to be a 
pseudo-random matrix of the same type since multiplying 
Fkr by itself M times is expected to establish correla- 
tions between the matrix elements (mi \F^ R \m2)- Nev- 
ertheless, we do assume that the matrix {m\FMKR\m') 
is banded and pseudo-random since the eigenstates of 
Fmkr have no simple connection with those of Fkr- 
Consider now an arbitrary matrix element (mi\F K I R \m2) ■ 
Due to the quantum diffusive dynamics within the M 
kicks, the (milF^Jm^) with \m\ — m%\ » 1 should 
be much greater than (mi |-FVii| m^) (nevertheless, both 
of them can be very small). According to Eq. JSJ, 
this implies that the matrix element (mi\FMKR\'m,2) is 
also much greater than (mi\FKR\m2)- In this sense, we 
expect that the {m\FMKR\m') matrix should display a 
wider band than does {m\FKR\m'}. 

The band width of the matrix (m\FKR\m') can be de- 
fined by choosing a cut-off value for its matrix elements. 
One traditional choice is |mi — m,i\ ~ k. In this case, 
(mi \FKR\m2} at the boundary of the band is on the order 
of Jkik), a number which is sufficiently small. However, 
this cut-off value, if applied to the matrix (m\FMKR\m'), 
would yield almost the same band width as Fkr- Hence, 
quantitatively characterizing the band structure of the 
matrix (m\F]\tKR\m') is subtle, since the very small ma- 
trix elements (mi \FMKR\m2} must play an important 
role in enhancing the dynamical localization length of 
MKR. 

To gain more insight we computationally examined the 
M-dependence of the dynamical localization length of the 
MKR. Specifically, we numerically diagonalized Fmkr, 
where each matrix is generated using 2 14 basis states and 
is then truncated at dimension d chosen below. We char- 
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FIG. 2: (a) The M-dependence of the dynamical localiza- 
tion length of MKR, characterized by the Shannon entropy 
Smkr averaged over all approximate eigenfunctions of the 
MKR propagator Fmkr- (b) The M-dependence of the band 
width bMKR associated with the MKR propagator Fmkr, 
defined by a cut-off value (as small as 10 -20 ) for its matrix 
elements (mi\FMKR \m2). 



acterize the average dynamical localization length by the 
Shannon entropy Smkr |27j| (note that the Shannon en- 
tropy is simply proportional to the dynamical localization 
length [23) averaged over all approximate eigenstates 
\4>i) of Fmkr, i.e., 



S 



MKR 



J^exp 



d/2 



|H0,)| 2 ln|H<fe-)| 2 



=-d/2 



where the constant a equals 0.96, and d is chosen to 
be 2700. The results for the M-dependence of Smkr, 
for M = 2 to M = 400 is shown in a log — log plot 
in Fig. Smkr is seen to behave initially as a 

smooth increasing function of M, and then to saturate at 
M ~ 50. To explain the results from the perspective of 
Band-Random-Matrix theory, we choose a cut-off value 
for matrix elements (mi\FMKii\n r i2) so as to define the 
band width bMKR- Interestingly, we find that this cut- 
off value must be extremely small (roughly speaking, at 
least as small as 10 -10 ) in order that the M-dependence 
of bMKR resembles that of Smkr- For example, Fig. [2b 



displays bMKR as a function of M for a cut-off value of 
10~ 20 . The evident similarities between Fig. and Fig. 
|2t> suggest (i) that we can indeed qualitatively explain 
the enhanced dynamical localization length in terms of 
the Band-Random-Matrix theory, and (ii) that even ex- 
tremely small quantum fluctuations in the values of the 
matrix elements of the MKR quantum map operator af- 
fect its dynamical localization length. 

In particular, comparing Fig. [2^ with Fig. [2b, one 
sees that the saturation behavior of Smkr for large 
M is consistent with the saturation behavior of bMKR- 
The latter reflects the saturation of the matrix ele- 
ments (mi \Fmkr\tti2) and therefore the matrix elements 
{mi\F^ R \m2). As such, we infer that the saturation be- 
havior of Smkr is simply a result of DL in KR. 

Detailed studies on numerous other cases with varying 
k and r show that the above result, i.e., Smkr first in- 
creases with M and then saturates, is quite general, as 
long as the system has a completely chaotic classical limit 
and quantum correlations are insignificant. On the other 
hand, if Smkr behaves differently, then there are two 
possible origins: either the system is in the deep quan- 
tum regime or there are non-negligible regular islands in 
the classical phase space. For example, in the next sec- 
tion we show cases in which the energy absorption in the 
MKR with M = 2 is appreciably larger than that in the 
MKR with M = 3. In these cases one can obtain even 
more significant changes in DL. 



IV. CLASSICAL ANOMALOUS DIFFUSION 
ASSISTED CONTROL 

In the previous section we studied cases where both KR 
and MKR essentially have a fully chaotic classical limit. 
However, as shown below, the MKR can also display new 
non-chaotic classical phase space structures that are ab- 
sent in KR. In particular, regular islands with very inter- 
esting transporting properties can be induced in MKR. 
In such cases one can achieve even more dramatic alter- 
ation of the DL than that shown above. 

Consider first the classical MKR map Eq. © for 
M = 2. Interestingly, in this case there exists transport- 
ing trajectories that are different from those in KR. In 
particular, we have previously observed |2(i| that for n — 
(2Z2 + l)7r, trajectories emanating from L = (21i + 1)tt, 9 — 
±7r/2 will be shifted by a constant value [±(2/2 + l)7r] in 
L after each kick. This observation suggests that new 
transporting regular islands that differ from the acceler- 
ator modes in KR can be created by changing the sign of 
the kicking potential after every two kicks. This is con- 
firmed in our extensive numerical studies, both here and 
inRef. HJ. 

Note that in our previous work [20|. we were most in- 
terested in the quantum-classical comparison in anoma- 
lous diffusion and considered relatively small effective 
Planck constants. In that case we found that quantum 
anomalous diffusion induced by the new transporting reg- 
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ular islands can be much faster than the underlying clas- 
sical anomalous diffusion. Here, to make a closer con- 
nection to atom optics experiments we consider larger 
t ~ f .0. In these cases the effective Planck constant is 
about an order of magnitude larger than the area of the 
phase space structures associated with classical anoma- 
lous diffusion. Intuitively, one would anticipate that such 
transporting regular islands are too small to be relevant 
to the quantum dynamics. Surprisingly, this intuition is 
incorrect, as shown below. 
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FIG. 3: Classical phase space structures of (a) the standard 
map and (b) the map of Eq. Q with M = 2, in the case of 
k = 5.0. All variables are in dimensionless units. Note that 
the small regular islands seen in panel (b) are transporting 
while those in panel (a) are not. 

To be more specific, consider first the case of k = 5.0. 
Figure |21 displays the classical phase space structures of 
both KR and the MKR with M = 2. While the regular 
islands seen in Fig. |2K (the KR case) are not transporting 
(that is, the momentum of the trajectories launched from 
these islands is bounded and oscillates periodically), a 



simple computation reveals that the small islands seen in 
Fig. |2t> in the MKR case are transporting. That is, clas- 
sical trajectories launched from the right (left) transport- 
ing regular island shown in Fig. have their momentum 
shifted by tt (— 7r) on the average after each kick, indi- 
cating that these islands originate from the marginally 
stable point L = (2Zi + l)n,6 = ±7r/2 with k = ir. 
Hence, in this case phase manipulation in going from 
KR to the MKR has both destroyed the nontransport- 
ing regular islands of KR and induced new transporting 
regular islands. Consider a second case with k = 10.0. 
The corresponding classical phase space structures are 
shown in Fig. (KR) and Fig. [§) (MKR) (Note that, 
to clearly display the transporting regular islands, only a 
part of one phase space cell is shown here). While there 
are hardly any regular islands seen in Fig. two small 
transporting regular islands are seen in Fig. 0Jd. The av- 
erage momentum shift for each kick associated with these 
two islands is found to be ±37r, consistent with the fact 
that k = 10.0 is close to 37r. 

Consider now the quantum dynamics of these systems. 
There have been only a few studies on the quantum 
dynamics of delta-kicked systems where classical chaos 
coexists with transporting regular islands. Of particu- 
lar relevance is the previous result that the accelerator 
modes of KR enhance deviations from the normal DL be- 
havior in KR 28, 29], even for systems far from the semi- 
classical limit. Since MKR displays new transporting is- 
lands, we therefore anticipate that DL may be strongly 
affected by modifying the Hamiltonian from the KR to 
MKR system. This is indeed seen below. The results 
are, however, counter-intuitive, since the classical trans- 
porting regular islands created by phase manipulation of 
the kicking field are found to have an area that is much 
smaller than the effective Planck constant. 

For example, for each of KR and MKR, Figs. [5] and 
El display energy absorption for the cases of r = 1.0, 
k = 5.0 and r = 1.0, k = 10.0 (corresponding to K = 5 
and k — 10), respectively. Also shown is the MKR case 
with M = 3, discussed below. It is seen that the energy 
absorption associated with the MKR with M = 2 (upper 
dashed line) is much larger than that of KR (solid line). 
Consider, for example, £ at a specific time t = 1500T 
when the energy absorption of both KR and MKR has 
clearly shown signs of saturation (e.g., the average rota- 
tional energy may decrease with time due to statistical 
fluctuations). In the first case (Fig. [5J, E = 67.7 for KR 
and E = 1269.7 for MKR. In the second case (Fig. EJ, 
E = 798.0 for KR and E = 10049.0 for MKR. In both 
cases a control factor larger than an order of magnitude 
has been achieved in going from the KR Hamiltonian to 
the MKR case with M = 2. 

This is not the case for M — 3, r = 1.0, k = 5.0 shown 
in Fig. |SJ Here the energy absorption in the MKR with 
M = 3 is only slightly larger than in the KR and far less 
than in the MKR with M = 2. Similarly, for the case 
of t = 1.0, k = 10.0 shown in Fig. El although energy 
absorption in the MKR with M = 3 is much enhanced 
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FIG. 4: Same as in Fig. ^except k = 10.0. Note that the 
small regular islands seen in panel (b) are transporting. 



(compared with KR), it is still not as significant as in 
the MKR with M = 2. This indicates, as confirmed by 
directly examining P(m) (not shown) after saturation, 
that for both cases the dynamical localization length of 
the MKR with M = 3 is no larger than that of the MKR 
with M = 2, contrary to what is observed in the previous 
section. This is because the underlying classical dynam- 
ics of MKR here is not completely chaotic, i.e., the clas- 
sical dynamics displays characteristics of anomalous dif- 
fusion due to transporting regular islands in phase space. 
Thus one can expect statistical deviations from the Band 
Random Matrix theory, used previously to relate M to 
the extent of control. 

These results emphasize that the control mechanism 
here is uniquely based upon the transporting regular is- 
lands created by our control scenario. This is further 
supported by the lineshape for DL, which can be strongly 
nonexponential, as discussed below. 



FIG. 6: Same as in Fig. ^except r = 1.0, k = 10.0. 



V. NONEXPONENTIAL LINESHAPES FOR 
DYNAMICAL LOCALIZATION 

It was pointed out more than two decades ago p| that 
the DL of KR can be mapped onto the problem of An- 
derson localization in disordered systems. In particular, 
an exactly soluble case of disorder in tight-binding mod- 
els, i.e., the Lloyd model [3(j, suggests that dynamical 
localization should assume an exponential line shape, at 
least on the average. This has been confirmed by numer- 
ous computational studies on KR. For example, Fig. 
clearly demonstrates, for both KR and MKR, that the 
distribution function P(m) can be fit beautifully with 
an exponential function with a characteristic localization 
length. 

However, the Hamiltonian nature of the KR and MKR 
implies that there always exist some subtle quantum 
phase correlations in the quantum dynamics. Hence, in 
addition to some universal properties of DL, the DL line- 
shape can display rich non-universal properties, e.g., the 



7 



system may display nonexponential dynamical localiza- 
tion. Noncxponcntial dynamical localization has been 
previously observed in KR but its origins are still poorly 
understood [2ll3l|. 

Here, we demonstrate that the MKR with M = 2 
can display strongly nonexponential line shapes for DL, 
rarely seen in KR. We focus on the MKR with M = 2 
since the classical MKR with M = 2 has transporting 
regular islands that are absent in KR. Second, transport- 
ing regular islands may induce large fluctuations in DL 
|28|. However, we examine below cases with connections 
to anomalous diffusion as well as those without clear con- 
nections to anomalous diffusion. We have also studied 
other versions of MKR with and have found that 

nonexponential lineshapes for DL in the latter case are 
much less common than in the M = 2 case. 
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FIG. 7: Four examples of nonexponential lineshapes for dy- 
namical localization in the MKR with M = 2, shown in terms 
of the probability P(m) of finding the system in the state \m) 
after 8000 kicks, with the initial state |0). In each case the 
broad lineshape is associated with the MKR, and for the pur- 
pose of comparison, the narrow lineshape of the analogous 
KR is also shown, (a) r = 1.0, k = 5.0, (b) r = 2.0, k = 5.0, 
(c) t = 1.0, k = 5.7, and (d) r = 2.0, k = 6.0. 

Figure H compares nonexponential lineshapes for DL 
in the MKR (upper dashed lines) to the analogous expo- 
nential lineshapes for DL in KR (solid lines), for four dif- 
ferent values of k and r. The lineshapes are obtained by 
propagating the quantum dynamics for 8000 kicks from 
the initial state |0), and will remain essentially the same 
for longer propagation times 32]. The huge difference be- 
tween the KR and MKR lineshapes is striking. As shown 
in Fig. H in the MKR case, P(m) plotted in the log scale 
displays structures that are far from a purely exponential 
lineshape. For example, one sees that the initial exponen- 
tial decay rate of P{m) with \m\ is considerably smaller 
than its large- \m\ exponential decay rate, suggesting that 
multiple characteristic lengths are needed to describe the 
DL of the MKR. Also seen is that the difference in P(m) 



between KR and MKR can be as large as ten orders of 
magnitude or more. Further, Figs. H> and [7]: indicate 
that the decay rate of P(m) for large |m| in the MKR 
case is very similar to the KR case, although this is not 
the case in Figs. Hi and 01. 

To the best of our knowledge, this is the first demon- 
stration that intriguing differences in the lineshapes for 
DL can be created by simply changing the sign of the 
kicking potential periodically. 

The four cases of nonexponential lineshapes shown in 
Fig. [7|can be divided into two categories, based upon the 
properties of their underlying classical dynamics. The 
classical dynamics associated with the cases in Figs. [7^ 
and 03, shown in Figs. |2t> and 03, displays transporting 
regular islands. The presence of these classical structures 
implies an inhomogeneous classical phase space, and, as 
demonstrated in the previous section, may have a signif- 
icant impact on the quantum dynamics even when their 
size is much smaller than the effective Planck constant. 
In this regard, the nonexponential lineshapes shown in 
Figs. 0, an d 03 may not be totally surprising. How- 
ever, for the other two MKR cases shown in Figs. and 
Hi, we did not find any regular islands in their classical 
phase space even when examined on a very fine scale, sug- 
gesting that their classical dynamics is essentially fully 
chaotic. Thus, understanding the nonexponential line- 
shapes shown in Fig. and Fig. Hi will be even more 
challenging. 

The nonexponential lineshapes for DL arise from ex- 
tremely small quantum fluctuations and residual quan- 
tum correlations in quantum chaos. To be able to resolve 
the nonexponential lineshape for DL, P(m) has to be 
known with high precision. For example, the two shoul- 
ders shown in Fig. 03 involve occupation probabilities 
P(m) as small as ~ 10 -14 . It is therefore not surprising 
that, while it is common to have strongly nonexponential 
lineshapes for DL in the MKR case, each individual line- 
shape is highly sensitive to the exact value of the effective 
Planck constant. For example, for the case shown in Fig. 
03, increasing the value of r from 2.0 to 2.0 + 10~ 5 can 
completely destroy the nonexponential lineshape! This 
drastic change in the lineshape for DL even occurs with- 
out causing an obvious difference in energy absorption 
behavior. Evidently then, both experimental observa- 
tions and theoretical predictions of nonexponential dy- 
namical localization are far from trivial and are in need 
of further study. 



VI. DISCUSSION AND CONCLUSIONS 

This paper has dealt with control of dynamical local- 
ization in kicked rotor systems. In all cases we manipu- 
late the external kicking field to alter the properties of the 
rotor system, i.e. the distribution of population amongst 
rotor energy levels after saturation as well as the energy 
absorption. In particular, we have examined the effect of 
introducing a reversal of the kicking field after M kicks 
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which, within the framework of quantum mechanics, cor- 
responds to introducing a phase shift amongst rotor en- 
ergy levels. 

Two parameter regimes have been examined, one 
which shows enhanced DL lengths with increasing M, 
and the other which need not. This behavioral differ- 
ence can be understood in terms of the character of the 
underlying classical phase space: the former systems are 
completely chaotic whereas the latter show a mixed phase 
space that includes transporting regular islands. Indeed, 
we have found that even if the transporting islands are 
tiny compared to the effective Planck constant, they still 
have a profound effect on the control of the DL. Fur- 
ther, a comparison of the traditional kicked rotor system 
to the modified kicked rotor systems shows that the lat- 
ter is much more capable of displaying nonexponential 
dynamical localization. Thus, by modifying the kicking 
potential are able to control the dynamical localization 
in the kicked rotor. 

The results of this paper are relevant to two fields 
of study: quantum control and kicked rotor dynamics. 
From the control perspective, modifying the kicking field 
changes the dynamics. However, this system does not 
obviously permit a picture in terms of interfering quan- 
tum pathways (the standard view of weak field coherent 
control [lR fl^|) since (a) the kicking field is always on, 
and (b) there are a multitude of interfering transitions 
responsible for the observed behavior. Indeed, it is even 
difficult to isolate the interfering pathways that are re- 
sponsible for dynamical localization in the simple kicked 
rotor, whose dynamics is easier than that of the MKR. 

From the viewpoint of kicked rotor studies in the 
field of quantum chaos, this paper provides insights into 
the quantum dynamics in the case displaying classical 
anomalous diffusion. As one of the results of this study, 
we find that classical transporting regular islands can 
dramatically affect the quantum dynamics even when 
their size is much smaller than the effective Planck con- 
stant. Further, the MKR system proposed in this pa- 
per provides a new model for the study of quantum dy- 
namics where the underlying classical chaos coexists with 
transporting regular islands. By choosing proper system 
parameters, we can create transporting regular islands 
whose size varies from being much larger to being much 
smaller than that of the accelerator modes of KR. En- 
couraged by this, we plan in the near future to further 
use the MKR to study quantum tunneling between the 
transporting regular island and the chaotic sea [2(1 Eflj ] 
and between transporting regular islands, and to study 
the phase space structure of quantum eigenstates |24| . 

In the fully chaotic case, our approach suggests the 
need for additional studies of dynamical localization from 



the Band-Random-Matrix theory perspective. For exam- 
ple, we have qualitatively explained the results in Sec. 11111 
in terms of well known features of Band- Random-Matrix 
theory. However, we found that the band width of the 
quantum map operator had to be defined using an ex- 
tremely small cut-off value for the matrix elements, sug- 
gesting that a quantitative understanding of the MKR 
results such as Fig. [21 may require new models of Band- 
Random-Matrix ensembles. 

The strongly nonexponential lineshapes for DL found 
in the MKR with M = 2 further demonstrate the need 
for more theoretical work on properties of DL. In partic- 
ular, our results should motivate greater interest in char- 
acterizing and understanding nonexponential dynamical 
localization, with efforts directed towards explaining why 
nonexponential dynamical localization occurs for some 
system parameters and not for others. This is of im- 
portance in understanding the high sensitivity of non- 
exponential lineshapes for DL to the exact value of the 
effective Planck constant. 

We have chosen the system parameters to be within 
the reach of current atom optics experiments on KR 0] . 
Although experimental studies of nonexponential line- 
shapes for DL are difficult, we believe that it is straight- 
forward to experimentally observe the results of Sec. IIIII 
and Sec. IIVI Apart from the atom optics realization of 
KR and MKR, it is also interesting to consider a molec- 
ular version of KR and MKR, i.e., diatomics periodically 
kicked in strong microwave fields 0, lisf . Preliminary 
computational studies |3^] confirm that directly observ- 
ing quantum control of dynamical localization in molec- 
ular rotational motion is possible, e.g., in the case of 
quantum anomalous diffusion. Another promising ex- 
perimental realization of KR and MKR requires kicked 
particles in a square- well potential [3^| . Along this direc- 
tion an interesting model (which is very different from 
ours) has recently been proposed for the study of classi- 
cal and quantum anomalous diffusion [35l |. 

In summary, consideration of control in classically 
chaotic quantum systems is of general interest and im- 
portance to both fields of quantum chaos and quantum 
control. In this paper we have demonstrated, via a modi- 
fied kicked rotor model, that dynamical localization, per- 
haps the best known phenomenon in quantum chaos, can 
be modified over a wide range. The results are of both 
experimental and theoretical interest. 
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